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1. Introduction 

In Bose-Einstein condensation, atoms with integer spin ('Bosons') occupy a single- 
particle quantum state macroscopically. This phenomenon, predicted by Bose [1] and 
Einstein [2] for non-interacting particles over 80 years ago, was observed experimentally 
in trapped dilute atomic gases in 1995 0,H]. Many similar experiments have followed [5]. 
These observations have renewed theoretical interest in the Bose-Einstein condensation 
of systems that lack translational invariance. For recent reviews, see e.g. [SI [3, [8]. 

Theoretical studies of Bose-Einstein condensation for zero temperature often, 
though by no means always j9j QUI [TT], make use of a macroscopic wavefunction $(x, t) 
that satisfies a cubic nonlinear Schrodinger equation [91 [10l [121 03] • Most recently, a 
mathematically rigorous derivation of this equation was given in the limit of an infinitely 
large number of interacting particles [Hj. 

The use of the nonlinear Schrodinger equation for $(x, t) has been deemed ade- 
quate for many experimental situations at extremely low temperatures [5]. However, 
this description is fundamentally incomplete even at zero temperature for a finite num- 
ber of interacting particles. In principle, atoms are excited from $ to other states. This 
many-body process is not accounted for by the nonlinear Schrodinger equation. 

Particle excitations in Bose-Einstein condensation were described systematically 
by Lee, Huang and Yang [15J for systems with translational invariance and periodic 
boundary conditions. In this setting, where there is no external potential, atoms are 
primarily excited in pairs from the (lowest) state of zero momentum to states with 
opposite momenta [15] . This process leads to phonons and sound vibrations [T5] . 

Wu [9] extended the theory of [15J to systems that lack translational symmetry. 
A key ingredient of his formulation is the pair-excitation function, Kq(x, y, t), which 
describes the scattering of atoms in pairs from $ to other states at positions x and 
y; see (|2]) below. This formulation yields coupled nonlocal equations for $(x, t) and 
fT (x, y,t) [9j [10], and therefore transcends other treatments based solely on the non- 
linear Schrodinger equation. However, solutions of the resulting coupled equations have 
remained largely elusive. By approximately decoupling the two equations, Wu described 
a time-independent solution for K Q in a slowly varying trapping potential [10] . 

In this paper we address analytically aspects of the time dependence of Kq for a 
system of trapped interacting atoms at zero temperature. Our starting point is the 
coupled system of $ and Kq [HI [TO] , which we simplify by assuming that $ satisfies the 
nonlinear Schrodinger equation independently. The resulting nonlocal equation for K 
is solved approximately for a time- independent, slowly varying external potential. In 
this context, $ is treated as a given, variable coefficient of the equation for Kq. We 
carry out a large-t asymptotic analysis that shows how Kq approaches the steady state 
if Kq = at t = 0, i.e., initially all atoms occupy the single-particle state 

Our study of pair excitation is motivated by three broader questions. The first 
question concerns the precise manipulation of atomic gases at very low temperatures. 
In current experimental setups, a rich variety of effects are observed including the de- 
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pletion of the macroscopic state <3>; see e.g. [TB]. Therefore, it is of interest to refine our 
understanding of how atomic interactions influence macroscopic properties of the Boson 
gas. The second question concerns the possible extension of the concept of a phonon 
to systems that lack translational symmetry. The third question concerns the analysis 
of coupled nonlocal equations. Pair excitation leads to integrodifferential equations of 
motion that are as yet unexplored. Here, we seek particular solutions but abandon 
mathematical rigor. Analytical solutions of simplified models such as the one in this 
paper may be a guide for future rigorous studies of similar models. 

The starting point is Wu's formulation [9j [10] for a system of N pairwise interacting 
Bosons at positions {x,}. By units with h = 2m = 1 (h: Planck's constant, m: atomic 
mass), the many-body Hamiltonian reads 

N d 
H = Y^[-&i+Ve(xi)] +4na'%26(x i -x j )—x ij , a>0, (1) 

i=l i^j 3 

where Aj is the Laplacian corresponding to Xj, a is the scattering length, V e is the 
external potential, and := |xj — By use of quantized fields [IT], the iV-body 
wavefunction of the Boson system is assumed to be of the form [9j [10] 

q?(t) = Af(t) e v{t) {(N\y 1/2 a*(t) N \vac}} . (2) 

In this equation, M is a normalization constant, a^t) is the creation operator for the 
state | vac) is the vacuum state, and the operator V(t) is 

V{t) ex j J dxdy^(x,t)^(y,t)^o(x,y,t)a (t) 2 . (3) 

Here, ip*{"X-, t) is the creation field operator corresponding to the space orthogonal to 
$. Notice that the integrand of V in ([3]) describes the annihilation of two particles 
from $ with the simultaneous creation of two particles at positions x and y at other 
states. If V(t) = (or, K = 0) then no atoms are excited from $. So, the many-body 
wavefunction reduces to a tensor product of single-particle states each of which is 

The functions $(x, t) and K (x.,y,t) are found to satisfy two coupled nonlinear 
equations [9l[T0]. These equations take the form [9[ [TO] 

i dt$ = [~A X + 7 e (x) + 8vrap |$r - 4vrap C(t)]$ + iV- x Wi(x, t), (4) 

[idt - 2E(t))K = -(A x + A y )K + Svrapo $(x, t) 2 5(x - y) 
+ {-Z(t) + K(x) + V e (y) + 16vrap [|$(x, t)\ 2 + |$(y, t)\ 2 ]}K 

+ 8vrap / dz $*(z, t) 2 K (x, z, t)K (y, z, t) + iV- 1 W 2 (x, y, t), (5) 

where p is the (constant) equilibrium density of the system and i 2 = — 1. The functions 
C(t), Z(t) and E(t) are integrals of $ to be defined in section [2j and W/ (/ = 1, 2) are 
nonlinear functionals of $ and Kq. We note in passing that po can be eliminated from 
the equations of motion by using the variable $ := p $ in place of $ [JO]. The system 
of (jl]) and dSJ) form a nontrivial extension of the cubic nonlinear Schrodinger equation 
for $, which results from (pO) by setting Wi = 0. 
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In this paper we study (J5J) when V e (x) is sufficiently slowly varying by enforcing 
Wi = W2 = and $(x, t) = e~ tEt Q(x.), where E is the energy per particle of the 
macroscopic state. Despite these crucial simplifications, (jSJ) remains nonlocal and cou- 
pled to the solution of the nonlinear Schrodinger equation. Equation (jSJ) can be solved 
exactly in terms of a spherically symmetric Fourier integral. We evaluate this integral 
for sufficiently large t for any value of r = |x — y| by expanding the integrand in an 
appropriate series and applying the standard stationary-phase method [IS] . As t and r 
vary, significant contributions to Kq occur when a stationary-phase point coalesces with 
the origin in Fourier space, which is the endpoint of integration. It is of some interest 
to recognize this coalescence as an 'elementary catastrophe', by analogy with studies 
of diffraction phenomena by Berry [19]. The time-dependent pair-excitation function is 
described by a solution of Lommel's differential equation [20J. 

We stress that our work relies directly on ansatz (j2J) for the many-body wavefunc- 
tion: the pair-excitation function, if (x, y, i), is used as a variable in addition to the 
macroscopic wavefunction, $(x, t). This approach is distinctly different from methods 
based on the hydrodynamic theory of superfluids, e.g. [2D [22] . Once K and $ are 
known, many-body properties of the Boson gas can in principle be computed via ([2]). 

A few remarks on our main simplifying assumptions are in order, (i) K Q does not 
act back on $ since we take Wi = in (J4j). Thus, the issue of how $ is modified by 
Kq is left unresolved here. Decoupling $ from Kq in this sense may pose a limitation 
on the time scale for the validity of our results in a physical setting. Therefore, the 
large-t limit studied here should be interpreted with caution, especially in connection 
to experiments, (ii) The external potential, V e , is time independent and slowly varying, 
and $ is taken to be time independent. A reasonable approximation for $ then results 
by neglecting the Laplacian in the nonlinear Schrodinger equation [10J. This simplifi- 
cation, sometimes referred to as the 'Thomas- Fermi approximation' in the context of 
Bose-Einstein condensation [6j [7] , amounts to seeking an outer solution in the sense of 
singular perturbation [TSJ. The effect on the motion of K of possible 'boundary layers' 
for $, where the Laplacian needs to be retained pi)], is not addressed by our analysis. 

The remainder of the paper is organized as follows. In section [2] we review the 
simplified equations of motion. In section [3] we focus on the instructive case with 
translational invariance where V e and $ are constants: in section [37X1 we derive an integral 
formula for K (x, y, t); and in section [3721 we evaluate the requisite integral for large t. 
In section H] we focus on the case with a slowly varying external potential: in section I4TTI 
we describe an approximate solution for $(x, t) revisiting [ID]; and in section I4T21 we 
show that the time-dependent Kq can be determined with minor modifications of the 
analysis of section [31 In section [5] we summarize our results and discuss related open 
questions. The units with h = 2m = 1 are used throughout the analysis. 
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2. Background theory: simplified equations of motion 

In this section we review the simplified equations of motion for $ and Kq. which are 
derived in [10] . These equations form the starting point of the analysis in sections [3] 
andU 

The macroscopic wave function $(x, t) satisfies the nonlinear Schrodinger equation 
tdM^t) = [-A x + K(x) + 87rap |$(x,t)| 2 -47rapoC(t)]$(r,t) , (6) 
where po = N/Q, Q is the volume of the system, 

C(t):=n -1 / dx|$(x,t)| 4 , (7) 
and $ is subject to the normalization condition 

n~ x J dx|$(x,t)| 2 = i . (8) 

The pair-excitation function K (x., y, t) solves the integrodifferential equation [10] 
[id t - 2E(t)]K (x, y, t) = -(A x + A y )K (x, y, t) + 8vrap $(x, tf 
x 5(x - y) + {-2£(t) - 16vrap ((t) - 2( e (t) + V e (x) + V e (y) 

+ lQnap [|$(x,t)| 2 + |$(y,t)| 2 ] W (x,y,t) + 8vrap 



x J dz$*(z,t) 2 K (x,z,t)ir (y,z,t) -87rapo^ 1 {$(x,t)$(y,t) 
x [|$(x,t)| 2 + |$(y,t)| 2 -C(t)] + <f(x,t) J dzK (y,z,t)|$(z,t)| 2 
x$*(z,t) + $(y,t) j dz^o(x,z,t)|$(z,t)| 2 $*(z,t)} , (9) 

where 

E(t) :=n~ l J dx$*(x,t) id t $(x,t) , (10) 
C(t):=Cr x J dx|V$(x,t)| 2 , (11) 

Ce(t) J dxK(x)|$(x,t)| 2 . (12) 

For sufficiently large N and fixed po, (EJ) is further simplified: the term proportional 
to Q _1 in the right-hand side is neglected. Thus, (EJ) reduces to the equation [TO] 

[i d t - 2E(t)]K (x, y, t) = -(A x + A v )K (x, y, t) + 8vrap <l>(x, tf 

x 5(x - y) + {-2C(t) - 167rap C(t) - 2C e (t) + K(x) + K(y) 

+ 167rap [|$(x, t) | 2 + |$(y, t) | 2 ] Wx, y, t) 



+ 87rap / dz$*(z,t) 2 ir (x,z,t)ir (y,z,t) • (13) 
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Following [TU] we use the center-of-mass coordinates, 




(14) 



and set 



#o(x,y,t) =: £(R,r,*) . 



(15) 



Hence, ( ITBl) is recast to the equation 



[i d t - 2E(t)])C(R, r, t) = -(jAfl + 2A r )/C + 87rap $(R, t) 2 5(r) 
+ { - 2f(t) - 167rapoC(t) - Ke(t) + V e (R + |r) + 7 6 (R - §r) 
+ 167rapo[|$(R+|r,t)| 2 + |$(R- |r,t)| 2 ]}/C(R,r,t) 

+ 87rap / dw $*(R - |r - w, t) 2 /C(R + |r - |w, w, t) 



x /C(R — |w, w — r, t) . 



(16) 



In the next section, we focus on solving (jSJ) and (|T6|) . 
3. Translational invar iance 

In this section we enforce translational invariance and periodic boundary conditions 
taking V e = const. First, we set $ = const, and study the resulting equation of motion 
for Kq. By invoking the Fourier transform of K (x., y, t) in x — y we find an integral 
representation for Kq. Second, we evaluate this integral for sufficiently large t. The 
results obtained here are useful in section H] where V e {x) is taken to be slowly varying. 

3.1. Solution for pair excitation 
In view of condition (jSJ) we set 



/C(R, r, t) should be R-independent for t > 0. With the replacement /C(R, r, t) = /C(r, t), 
(JIB]) becomes 

i d t )C(r, t) = -2A r /C(r, t) + 8vrap 6{r) + 167rap /C(r, t) 



$(x,t) = l t>0. 
By flZD, (HO]) and (JUJ), we have 

COO = 1 , E(t) = , ((t) = . 
Thus, (EJ) along with (TT2"j) entail 

V; = -47rap = Ce(*) • 

The goal is to solve (03} under ([T7j)-(fll)]). 
By the initial condition 

£(R,r,0) = /(r) , 



(17) 



(18) 



(19) 



(20) 




(21) 
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In addition, we have the symmetry relation fC(r, t) = fC(—r,t) assuming that /(r) = 
/(— r). Thus, a convolution integral appears on the right-hand side of (j2ip . 
Accordingly, we introduce the Fourier transform of /C by 

£(k,*) = J dr /C(r,t)e" ik - r , (22) 

whence 

f Hk ^ 

/C(r,t) = y ^e*-^(k,t) . (23) 

Hence, (1211) is transformed to 

i d t £(k, t) = SirapvK 2 + 2(A; 2 + 8ixap )K, + 87rap (A; = |k|) . (24) 

The steady state for ( 12T|) stems from taking = in ( |24"1) . We interpret the 
resulting time independence as the long-time (t — > +oo) limit, in an appropriate sense, 
of the solution to ( 1241) . The corresponding limit for JC reads 



P ( k 2 + 87rap -k^k 2 + 16vrap ~ , n , oc .x 
iqk, +oo) = — w - =: g Q (k). (25) 

87rapo 

This time-independent solution is derived in [TU] for a slowly-varying external potential. 
The inversion of this ~go(k) gives [10] 

g {r) = n- 2 (4na Po ) 3/2 x(r)- 1 lm[S 0A (i X (r)) - S Oi0 (i X (r))}, (26) 

where go{r) := /C(r, +oo), r = |r|, S 1 ^ denotes the Lommel function [23J, and 

X(r):=(167rap ) 1/2 r. (27) 

We turn our attention to the time-dependent /C. Equation (1241 can be solved 



explicitly. Details of the integration procedure are given in Appendix A In view of (|20p . 
the solution reads 

£ ( k, t) = U k) - [i - urn 1 _ ^;~^L m , , <*> 

where 

:= fc^fc 2 + 167rap = [?o(A0 - ?o(^)~ 1 ] , (29) 

p(k) := ^)-f(k) , (30) 
l-?o(^(k) 

and /(k) is the Fourier transform of the initial data, equation (1201) . Note that 
— 1 <^o<OifA:>0. The function u(k) gives the phonon spectrum, or the energy of 
excitation from zero momentum into a state of momentum k [T5] . 

We next apply spherically symmetric initial data, /(r) = fo(r) in (|20I) . Thus, p(k) 
is replaced by the function 

p(k) = p (/c) := ^777777; > ( 31 ) 
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where fo(k) is the Fourier transform of fo(r). Equation (I28p becomes 

^o^w-n-aw ^.^^, . (32) 

We note in passing that the right-hand side of (132!) can be continued analytically in 
the complex fc-plane where singularities are located only in the first and third quadrants. 
We refer the reader to Appendix B for details when / = 0. This observation is relevant 
to the asymptotic calculation of section [3721 and the long-time limit fl25|) : by deforming 
the inversion path for K, slightly below the positive real fc-axis, we can expicitly ensure 
that e~ 2laJ( - fc )' — > as t — > +oo. Consistent with definition (1221 . this limit can be 
generalized to include solution (128]) for a reasonably wide class of initial data /(k). The 
limit (125]) is then readily recovered. 

The inversion of /C(k, t) is expressed by the integral 



i-oo 

JC(r, t) — — — / dk ksm(kr) 
2n 2 r J 



x \go(k) - [1 - go(kf] : • (33) 



l-g (k)p (k)e~ 2i ^ k ) t 
We have not been able to compute this integral exactly in simple closed form in terms 
of known special functions. 

Equation (1331) is written in the convenient form 

A(r,f) := (16nap )-^[}C(v,t) - /C(r, +oc)] = _ ^Ef!!l 

2ir z r 

+oo 

p (k) e" 2 *"^)* 



x y ^^wii-awi ^r- aMt)t . (34) 



where A is non-dimensional and /C(r, +oo) = go( r ) is defined by ( |26i) . The function 
A(r, t) expresses the deviation of the pair excitation function, /C, from the steady-state 
solution, go(r). 

In the following, we focus on zero initial data, i.e. /o( r ) = 0, by which all particles 
initially occupy the state <3>; cf (T2j). Hence, we have 

p (k) = g (k) . (35) 

By the change of variable k = (167rapo) 1//2 sinh?7, ( !M|) becomes 

+oo 

\ r g-4r? e -irsinh(2?7) 

A M } = J dr l [smh(2r / )] 2 sin(fsinhry) — — i — (36) 



where the non-dimensional variables r and r are defined by 

f := (167rapo) 1/ ' 2 r , r '■= {I6irap )t . (37) 



As discussed in |Appendix B[ the path of integration in ( l36l) can be deformed slightly 
below the positive real axis since there are no singularities in that part of the 77-plane. 
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3.2. Asymptotic analysis for t ^ (Wnapo) 1 

In this subsection, we determine A(r, t) for zero initial data and any r when t ^> 
(167rapo) _1 - In view of (13"4"1) . a key idea is that, if |po(^)| < 1 ; we can expand 
[1 — 'go(k)p (k)e~ t2u} ^ t }~ 1 in a geometric series. Thus, we have 

+oo 

£(M) -&(*) = (47rap )- 1 a;(A ; )^[? Q (A ; )po(A ; )]'e- ii2 -W t , (38) 

where we used ( )29|) . When both t and r are large enough, each corresponding integral 
under the summation sign for A(r, t) can be evaluated by the method of stationary 
phase |18j . As t and r vary, a critical point coalesces with the endpoint of integration, 
k = 0, giving rise to a behavior described by the one- variable Lommel function [23J. 

We proceed to carry out this program by resorting to (1361) . using the non- 
dimensional variables r and f where r 3> 1. The dependence on r, t and apo can 
be restored via definitions ( l3"7|) once the calculation has been completed. 

3.2.1. Case f = 0(1). When f is fixed, the major contribution to integration in (1361) 
comes from the vicinity of rj = 0. We obtain the expression 

+oo 



AM) ~ — ^ / d?7 {2r)Y{fr]) 



.2 /- ^ e 



-ir (27?) 



27r 2 f / v ; 1 - e~ iT ( 2r >) 







' 4^(1) = ^^, (39) 



8tt 2 r 4 r w 120 r 4 
where V^- 2 ) is the logarithmic derivative of the Gamma function and = 7r 4 /15 is 

the third derivative of ip{ z ) at z = 1 |24j . Formula (1391) is not particularly informative 
because of the absence of r-dependence. 

It of some interest to re-derive (1391) in light of (1381) . By writing 

M-l 



[I _ g-4r? e -jrsinh(2r?)j-l _ 



g— 4Zr? g— iir sinh(2»?) 



Z=0 



e -4Mr? g -iMtsinh(2»?) 



-[ g-4»? g-iTsinh(2r?) ' — ' (^0) 



A(r, i) from (1361) becomes 

1 



A(r,t) 



27T 2 f 



^/i(f,r) + fl w (f,r) 



z=i 



(41) 



where 



+oo 

7 z (f,r) = y drj [sinh(2r/)] 2 sin(r sinhr?) e" 4 '" e -^sinh(2r?) , (42) 
o 

+oo 

/ e -4(M+l)r? e -i(M+l)rsinh(2r?) 
dr? [sinh(2r ? )] 2 sin(f sinh^) - _ — e _ trsinh(2r?) . (43) 



o 
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By deforming the integration path into the lower 77-plane, we have 

D 2i(M+l) (2y) -(M+l)r (2y) 



r e 2i(M + L) (2y) e -(M + 

R M (f, r) ~ 4r J dy y 3 1 _ e -r{2 y) 



-^rib m ((M + l)(l-2i/T)) ~ — A - t ^— ^ , (44) 

4r 4 ~ J\ I J J 2r 4 (M+ l) 3 ' V 7 



uniformly in r, where M ^> 1. The integral l\ of ( )42l) is 



+00 

3f 



J,(f, r) ~ / dry (2r/) 2 (r~r/) e"^ = ^ , f = 0(1) , (45) 







for any I > 1. So, by (14T1) we have 

£71 + ^ 



. / x 1 3f 
AM) 



2vr 2 2r 4 



/ 4 
z=i 



^r" 4 , M^+oo. (46) 



In the following, we invoke ( j4T~i) and f j42l) when f is comparable to or larger than r. 

3.2.2. Case f > 0{It). In this case we need to describe the r-dependence of A by 
appropriately accounting for sin(fsinhr/) in (1361) . For this purpose, we apply f|4T!) for 
M -> +00: 

1 00 

A M) = ^5>( r ~' r )' ( 47 ) 

Z7r z r z — ' 
1=1 

where 

7,(f,r) =/ I)+ (f,r)- V(r,r) , (48) 



-00 



ij,±(f,r) := j. J drj [sinh(2r/)] 2 e" 4 ^ e -«™n*(2»/)±if smh„ _ ^ 


The integrals J^-t must be treated differently. Specifically, the major contribution 
to integration in Ii _ always arises from the neighborhood of the endpoint, 77 = 0: 

4 4 

Ij _(f, t) ~ — — ~ — all f > . (50) 

' v ; (21 t + r -iAlf (2Zr + f) 3 ~ v ; 

It follows that [25 

+00 1 

EU»~^/(l + ^ (51) 
1=1 

where ip"(z) := d 2 ?/>(;z)/d;? 2 . 

By contrast, computing I l)+ requires the evaluation of a stationary-phase contri- 
bution when f is comparable to r. Most of the remaining subsection is devoted to the 
analytical computation of Ii + . 
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Remarks on the integral To evaluate it is advisable to introduce the 

associated phase 

0(77) := (f — 21t cosh rj) sinh?y . (52) 
The stationary-phase points are the non-negative roots of the equation dQ/di] = 0: 

f cosh 77 — 21t cosh(2?y) = , (53) 
which is solved by rj = rji > where 



cosh^ > if A:=^>1- (54) 

Note that for n > 0, 

0(r#) > 0, 0"(^) := d 2 0/d??\ = (f - 8/ r cosh 77,) sinh^ < . (55) 

In particular, if 77/ > 0(1) then Q(rji) > 0{lr) and |0"(?7z)| > 0(It). In this case, the 
stationary-phase calculation is carried out as usual [18]; see ( l60l) below. 

However, caution should be exercised in the asymptotic evaluation of Ii t +(r, r) when 
the point rji = 77/ (f,r) lies too close to 77 = 0. This coalescence occurs when 77/ becomes 
of the order of |0"(^)| _1 ^ 2 while < Pi — 1 <C 1. By expanding rji near (3i = 1, we find 
that a more precise condition for this coalescence is 

\2lr-f\ = 0((/r) 1/3 ) . (56) 

Accordingly, we distinguish three main regions for (f, r) as shown below. In the end of 
this section, we exploit the overlapping asymptotic expansions for Ii in these regions in 
order to construct a connection formula. 

Region I: 2lr - f > (It) 1 / 3 . In this case, the stationary-phase point r\\ is not real. 
The major contribution to 7j + stems from 77 = 0: 



-OG 



7 li+ (f,T)~i I dr; (277) 2 e-^-^e- 4 ^ 





4 4 

(57) 



{2lr-f-iAlf (2/r-f) 3 

This result, combined with (148]) and (|50|) . furnishes (j4"5l) when f -C 2/r. 

We note in passing that if f < 0(2/r) yet 2/r — f ^> (Zr) 1 / 3 for all Z then we can 
sum up ( 1571) . Consequently, 

+ OO _ 

EU»-^f(l-^) • (58) 
1=1 

Formulas f[5"Tj) and 058)) combined yield 

1 1 



r \ , „ / r 



A M)~^~;a + • (59) 



2r/ V 2r/J 

This formula reduces to f)39p via the Taylor expansions of ± f/ (2r)) about 1. 
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Region II: f — 2lr (Zr) 1 / 3 . The standard method of stationary phase is now 
applicable. The integral Jj + is approximated by [18] 



r) ~ -^^e^ e i@ ^" [sinh(2^)] 2 , (60) 

where Q(r)i) and 0"(^O are given by (152]) and (155]) . 

It is of interest to take the limit of ( 160]) as A — > 1, i.e. ^ — > 0. With 

e( Vl ) ~ 2/r A (^ + - 2/r(^ + §77?) 

o5/2 

~^(A-l) 3/2 /r , (61) 

e"( Vl ) ~ -2 3 / 2 3 1 / 2 (A - l) 1 / 2 Zr as A - 1+ , (62) 
the stationary-phase formula (|60l) reduces to 



2 7/4 



2 5/2 



r) ~ ^£ (A - l) 3 / 4 exp [i^- (A - l) 3 / 2 Zr + . (63) 

Region III: \2lr — f\ = (^((Zt) 1 / 3 ). In this more demanding case, the stationary-phase 
point rji is too close to the endpoint of integration and Q(r]) must be expanded at rj = 0. 
Approximation (IBTl) indicates that terms 0(rf) must be retained in 

First, we consider the region 2h > f with 2/r-f = 0((It) 1 ^). With the definition 

J_ (2Zr - r~) 3 / 2 

7 '' 3 3 / 2 (/r)V2 lD4j 
and the change of variable r/ = 2(Zr) _1 / 3 (7^/2) 1//3 sinh(v/3) in (149]) . becomes 

+0O 



+ 0O 

7z y dt> [coshf — cosh(t>/3)] e 



o 

2 



3«/r 



o 



2 4z /2/r-f 



~\ 3/2 



—17; sinhi> 



3/ T + 3(— J 7 ' = 0(1) ' (65 » 

where is the Lommel function [23]. In the above, the e~ 4lT1 term has been neglected 
in the integrand. This simplification is adequate if, for instance, f = 0(t) and only 
values I = 0(1) are of interest. Note that if 2lr — > f for some Z then 7/^0 and the 
second term in fl65|) vanishes, as it should. 

Formula ( 165]) connects smoothly to approximation ( [57]) when 7; 3> 1. Indeed, by 
use of the asymptotic formula |23j 

1 8 

S 0> l{z) ~ ^-3 as |z| — > +00 , |Arg z\ < n , (66) 

the first term in (165]) is exactly canceled. Thus, (165]) reduces to (]57j) . 
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To account for large values of /, / = 0(r 1 / 2 ) by which 21 r — f = 0(1), the e~ 4lr> 
factor has to be retained in the integral (J4"9~j) for An inspection of fl65|) shows that 
7; has to be replaced by 

2 (2/r-f-i4/) 3 / 2 



: 3 3 / 2 (/r) 1 /2 
The ensuing approximation for J/ j+ in place of (1651) reads 

'2/r-f-i4/\ 3/2 



4 



2 4z 
3lr + ~3 



3/7 



V 



(67) 



(6* 



Next, we restrict attention to the region f > 2lr with f — 21 t = 0((l t) 1 ^ 3 ). This 
case is not essentially different from the previous one. Formula (1651) can be continued 
analytically to complex values of 7;, as 7z varies continuously from ji = | 7 ;| for f < 21t 
to ^le - * 371 "/ 2 or |7i|e l37r//2 when f > 2Zr. Both continuations should yield the same result 
for J; )+ since this integral is a single- valued function of 7;. Indeed, from Appendix C we 
have 



(-We-*"') [(•Tie-*/'')- 1 - S 0i| (i 7 ;e-^/ 2 )] 

= - i - i 7/ S„,i (71 = -* + ili S i (ji e°") 
= (7/e^ 2 ) [(irne"***)- 1 - S^t^e^ 2 )} . 



(69) 



Equation (1651) combined with (10. 5|) of Appendix C yields 



2 

3fr 



4 
3 



2W 



3/2 



3/r J 



SoAhtt + lVSe-^H^l) 



f >21t , 



(70) 



where, by (Jb 



2 • -- -i/ 3 < 
| 7 ,| =2-3- 3 / 2 (f-2/r) 3 / 2 (/r)- 1 /2. 
For |t{| 3> 1 and f — > 2Zr (i.e., $ — > 1), we show that (1701) connects smoothly to the 
stationary-phase contribution of ( 1631) . The leading term in the brackets of ( 1701) arises 
from i?{/3 according to the expansion [23] 



r(l), 
1/3 ^ 




7TX 



+ OO . 



(71) 



The next asymptotic term in (170|) comes from the S 1 and exactly cancels the —2/ (3/r) 
Thus, ( 1701) reduces to 



r, r ~ 



35/4 



r-2h 
It 



3/4 



3 t |7i|+ivr/4 



(72) 



This formula is in agreement with (163|) under the substitutions r = f3i (2lr) and 

| 7j | = 2 5/2 3 -3/2 (A _ 1)3 /2 /r „ T]^ — > . (73) 
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The inclusion of the e~ 4lr/ factor in (1491) for necessarily modifies (1701) . By defining 

u 2 (r~-2/r + ?4/) 3 / 2 
11 ' 3 3 / 2 (/r)V2 

~ |7/| [l + i(8/3)^ | 7 , 1" 1 ] 3 / 2 , 0<77,<1, (75) 



we replace ( I7D1 by the formula 



2 4 /f -2/r + i4/\ 3/2 

r <> + ~ "377 + 3 I 377 ) 

S^m + l^e^H^)] . (76) 

Connection formula for Next, we derive an asymptotic formula that connects 
the apparently disparate yet overlapping formulas ( !57l) . ( IBTil) and (l65l) for It is 

instructive, although not necessary, to distinguish the cases f < 2/r and f > 2/r. 
Assume that f < 2/r. By virtue of fl68l) . (r, r) reads 

. 2 4i (2W-r - i4l\ 3/2 

"3/7+3 I 37r J 

/2i (2/r-r~-z4/) 3 / 2 \ 

x ^ U (3/7p ) ■ (77) 

This formula can be continued analytically to the region f > 2/r as was shown above, 
but may break down when f — 2/r ^> 0((/r) x / 3 ) (in region II). Our task is to connect 
this formula with the stationary-phase result fl60|) by using a single expression. 

Let us now consider f > 2/r. By inspection of flTDl . f[T2"|) . fl?3"j) and fl75]) . we propose 
to start with the composite expression 

2 

37r 

where the constants C and C are to be determined. This formula must connect smoothly 
to the stationary-phase result (16"0"j) when r)i > 0(1) and, thus, Q(r)i) ^> 1. Recall that 
the phase function O(^) is defined by (|52|) . 

In order to find C and C we replace S^i in (ITSjl by its large-argument expansion, and 



~ -£ — - c s o4 (e(i»)[i + ^/s^e^rT/V-) , ( 7 * 



then match the asymptotic result with fl60|) . By |Appendix C[ we use the relation [23] 

S i{ze-™) = - S ,(z) - | v^e-/ 3 

— e lz+i7T/A - -, Id — 00, |Arg z\ < it . (79) 

2z z 

For rji > 0(1), is thus approximated by 




T ,~ s „ 2 C „ / 37T 



2 
3/7 



3/r 6(^) + Ul rji V 2[G(»7i) + i4Z tj z 

x gje(j7i)+i7r/4 e — 4i r?i 

3 /r 



2 efa 



C ^ / e i@ ^ )+-/ 4 e" 4 ^ . (80) 
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The term in the brackets of (I80p should be zero while the second term must be identified 
with (l60l) . Accordingly, we have the relations 

£ -4eM- < 81 > 



C=,/,|^[*W- (82) 



By use of the identity 

Q(r]i) f sinh^ — /r sinh(2?7i) (sinhr^) 



|0"(^)| 4/r sinh(2r]/) — fsinh^ 1 + 2 (cosh rji ) 



2 



we find 



1 sinhr?; , , , l0 . . 

^7fv/l + 2(co:h^ [sM ' (M) 

C=V3 - {cOShVl)3 . (85) 
y/1 + 2(cosh^) 2 

To further validate (1751) given (154"j) and (155|) . we consider f — 2/r = (^((/r) 1 / 3 ) (in 
region III) so that < rji < 1. In this limit, (El and (ESD entail C ~ (4/3)r? 3 ~ (2/3) M 
and C ~ 1, while the argument of 5 i in (1751) becomes approximately |7z|e~ J7r . Recall 
the definition ( 1731) of 7/. Thus, ( 1751) reduces to (1701) . 

We point out a refinement of our procedure. A comparison of (1751) with formula 
( 1761) reveals that the C of ( 154)) does not reproduce the expected prefactor of Lommel's 
function when f — 21 t = 0(1), i.e. 77/ = 0(r -1 / 2 ). A remedy to this discrepancy is to 
replace r\i in the C of (1541) by 



4i 

vr-=\lvi+^:- (86) 

This 7# stems from including the exponent — 4/r/ in definition (1521) for 0(r?). The equation 
for the resulting, modified stationary-phase point 17/ reads 

2/r cosh(2r/;) — fcoshrji — i4l = , (87) 

which has solution 

coshi)i = a + ygw.wr (88) 

Expanding this equation about Pi = I yields 

f{f A~l + »2/r l f + ^Z-2/r 

1 + + 3 = 1 + 3 2/r ' (89) 

by which tji ~ fji in view of ( 1561) . 

A connection formula for J; + follows from (175|) with (154)) - (1561) : 

2 (coshr^) 3 1 1 (sinhf//) [sinh(2f/z)] 2 



\/3 y/l + 2(cosh^) 2 Zr ^ ^1 + 2(cosh^) 2 
x S i (6(i»)[l +i(8/3)^e(^)- 1 ] 3 / 2 e--) . (90) 
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Asymptotic formula for A(r, t). A large-t asymptotic formula for A(r,t) 
(16%apo)- 3/2 lJC(v,t) - /C(r,+oo)] is obtained by flSD, (|48]), (Hp and (|90|): 



A(r,t) 



r/(2r)] 

E 



(cosh 77/ )" 



2vr 2 f [ ^ [^1 + 2(cosh^) 2 V3/r 



1 (sinhr^) [sinh(2^)] 2 
V3 v/l + 2(cosh^) 2 



+00 

+ E 

Z=[f/(2r)]+l 



2 

3lr 



4z (2lr - r -iAl\^ 12 



3h 



.2% (2lr ~r -iAlfl 2 

x 1 ¥ (377p 



-L<(l + f/(2r)) 



(91) 



where [x] denotes the largest integer that is less than or equal to x. The infinite series 
for A is absolutely convergent and can be simplified depending on the value of f/(2r). 
For instance, if r < 2r then ( 19TT) reduces approximately to the expression 

'2i (2r - r - 4i) 3 / 2 ~ 



AM) 



27T 2 r 1 3r 

1 

4^3 



4i^2r-f-4A 3/2 



2r 



3r 



r 

27 



3 

211- 



r 

27 



(3r)V2 

-3n 



(92) 



The last formula shows that the magnitude of A varies from 0(r~ 4 ) when f = 0(1) to 
0((fr) _1 ) as r approaches 2r. 

More generally, if 2(n — l)r < f < 2nr for some positive integer n = 0(1) then (!9~T!) 
is approximated by the following finite sums: 



AM) 



27r 2 r 



1 f n_1 

r 2f I 



1=1 



(cosh /ft) 



2 (cosh r]i) 2 \/3h 



1 (sinhr^) [sinh(2^)] 2 
v/l + 2(cosh^) 2 



+ 



x 5 0) i (efa)[l + 2(8/3)^ Gfa 



*-ll3/2 -itt 



+ 



3nr 
1 

4^3" 



4i / 2nr — r — Ai n 
3 \ 3nr 



3/2 



^ 1 + 



2r 



< 1 



r 

2~7 



2z (2nr - f - Ai n) 3/2 



(3nr)V2 



^E ' 



1=1 



r 

27 



(93) 



As r/(2r) approaches n the magnitude of A locally increases to become 0((nfr) _1 ). 

Recall that the dependence on r, t and apo is restored via definitions (I3TI) by which 
f = (le-Trapo) 1 ^ 27 " an d r = (16%apo)t. Note that the quantity (lGyrapo) 1 ^ 2 is the sound 
velocity in units where H — 1 — 2m [IS]. Equation ( 19TT) is further discussed in section [5j 
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The asymptotic analysis of this section can be extended to nonzero spherically 
symmetric initial data, /o( r ) 7^ 0, by which po(k) ^ %{k). The major contribution 
to the requisite integral ( 1381) still comes from the vicinity of k = 0. For example, for 
I = 0(1) we can apply the approximation \po(k) go{k)} 1 ~ [po(0) goiff)} 1 if Po(0) <?o(0) ^ 0. 



4. Slowly varying trapping potential 

In this section we use a potential V e (x) in (Q and ffl3l) that is slowly varying and in- 
creases with |x|. We invoke the approximate time-independent wavefunction $(x) for 
the macroscopic state derived in [ID] . Accordingly, we obtain the time-dependent pair- 
excitation function Kq by solving ( fl3i) . We show that, under certain conditions, the 
analysis for K is not essentially different from the procedure of section [3j 



4-1. Macroscopic wavefunction 

Here, we revisit briefly the solution of the nonlinear Schrodinger equation (J6]) that was 
studied in [TP] . 

By the replacement $(x, t) = e~ iEt Q(x.), the equation for $(x) is 

[-A x + K(x) + 8na Po $ (x) 2 - 4vrap C]$(x) = E $(x) , (94) 

where is the energy per particle of the macroscopic state and $(x) is chosen to be real. 
We consider a sufficiently slowly varying V e (x), formally described by V e (x) =: V e (ex) 
where < e < 1 and V e = 0(1). Accordingly, is neglected in ([91]). Thus, $ ~ $ 
satisfies the equation [10] 

[V e (x) + 8vrap $o(x) 2 - 4nap ( - £]$o(x) = . (95) 

The solution of this algebraic equation reads [10J 



. (87rap )- 1/2 [4vrapoC + S - ^(x)] 1 ^ , x € ft , 

$o(x) = < o , x i n , (96) 



where the region 1Z is defined by [10] 

ft = {x : 4vrapoC + £ > ^e(x)} . (97) 

The function $ must vanish as x approaches the boundary, <9ft, of ft from the interior 
of ft. Notice that $o by ( 1961) is also slowly varying. The multiplication of ( 1961) by $o( x ) 
and subsequent integration in x yield [10] 

£ = 47rap C + Ce • (98) 

The right-hand side of (1961) is not differentiable in x at dTZ. A remedy to this 
problem was provided in [10] by retainment of the Laplacian in (l94"j) . 
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4-2. Pair- excitation function 

In this subsection we determine the time-dependent pair-excitation function in the 
center-of mass coordinates under approximation ( 1961) . In view of integrodifferential 
equation (fT6l) . let 

/C(R, r, t) = e~ i2Et /C (R, r, t) . (99) 

In light of the translationally invariant case (section [3]), we treat R and r as slow 
and fast variables, respectively, in ( 1161) . This consideration amounts to neglecting the 
Laplacian Ar and eliminating the variable r in the arguments of V e and $ in f|T6|) . In 
addition, we have [10] 

/C (R + \r - |w, w)/C (R - |w, w - r) ~ /C (R, w)/C (R, w - r). (100) 
Thus, f lTBT) reduces to [TU] 

id t JC {R, r, t) 2A r /C + 8vrap $o(R) 2 ^(r) 

+ 2[-87rapoC - Ce + V e (R) + 167rap $o(R) 2 ]/C 

+ 87rap $o(R) 2 j dw /C (R, w, t) /C (R, w - r, t) . (101) 

The variable R enters (11 Oil) as a parameter. By virtue of (1961) and (1981) . the term 
in the brackets of (1 10 II) reads 

- 87rapoC - Ce + V e (R) + 167rap $o(R) 2 = 87rap $o(R) 2 • (102) 
Thus, (1101 1) results from (pT!) of the translationally invariant case with the replacement 
po p $o(R) 2 • (103) 

Therefore, the formulation and calculation of section [3] remain essentially the same here 
when $o(R) 7^ 0. The case with < l ) o(R) = leads to a simplified evaluation of K as 
shown below. 

In the following, we outline the main results in view of f l 1 03D . We solve (110 II) with 
the initial condition 

/C (R,r,t = 0) = /(r) , (104) 

where / is assumed symmetric, i.e. /(r) = /(— r). Thus, 

/C (R,r,t) = /C (R,-r,t) . (105) 

By f l25l) . fl28|) - fl30|) and (UOip . the Fourier transform in r of /Co is found to be 

£ Q (R, k, t) = §b(R, k) - [1 - ? (R, kf] 

p(R,k)e- 2 ^ R ' fc )* 
X 1-^ (R, k)p(R, k)e" 2 -( R ' fc )* ' (1 ° 6) 



where 



~ (Ti i\ 87rap $o(R) 2 

go(R,k):= , (107) 

P + 8vrap $o(R) 2 + fc v / ^+!HA(RF 
a;(R, k) = k^k 2 + 16vrapo<l»o(R) 2 , (108) 
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, , , g CR.k)-f(k) , s 

p(R,k) := — ; 109 

i-3b(R,*)/(k) 

The inversion of (11061) is carried out through integral (1231) . 

It is worthwhile noting that the /Co obtained from (11061) is not differentiable in 
R at the boundary dTZ. This feature is due to approximation (196]) . A remedy is to 
retain A# in the equation of motion for the pair-excitation function and describe the 
spatial changes of K Q along the local normal to dTZ. This analysis would also require the 
analogous modification of the approximation for $(R) described in [10J. These studies 
lie beyond the scope of this paper. 

With regard to R, we distinguish the following cases. 

Exterior of region TZ (R ^ TZ). Formula (1 1 6 [) is greatly simplified since $o = 0. We 
have 

A^(R,k,t) = f{k)e- 2ik2t , R^TZ. (110) 

Recall that TZ is defined by (1971) . By inverting ( 11 101) we obtain the R- independent 
function 

|r' — r] 



|2 



/exp I I 
dr' fir') v & 7 , R i TZ , (111) 



which is the solution of the linear diffusion equation in the three-dimensional space 
cf (UOip with $ = 0. This /C vanishes as t — > +oo for a reasonably wide class of initial 
data. In particular, if f(r) is of compact support with size L, t > L 2 and r > 0(\/i), 
then (111 II) yields the known similarity solution 

/C(R, r, t) ~ (Svrt)- 3 / 2 exp (-^) | dr' /(r') . (112) 
If / = then /Co by (1 11 II) vanishes identically outside 1Z. 

Interior of 1Z (R G 7Z). Equation (I106p is inverted through (12^1) but the integration 
result is not expressed in simple closed form. For spherically symmetric initial data, 
/(r) = f (r), /C (R,r,t) is given by (JM]) replacing /C by /C , co{k) by u(R, k), g (k) by 
p (R, k), and by 

/■d n ?o(R,fc) -/o(fc) 
Po(R,k) = ~ • (113) 

l-Sb(R,A0/o(*) 
For zero initial data, f = 0, we have po(R, k) = ^o(R> k). 

The asymptotic analysis of section 13.21 for fo = is applicable here as well. The 
non-dimensional variables (scaled distance and time) r and r are now defined by 

f(R) := [16vrapo$o(R) 2 ] 1/2r > r(R) := 16vrap <l>o(R) 2 t . (114) 

These definitions lead again to the integral (1361) where A is now defined by 

A(R,r,t) := (16vrapo$o(R) 2 ) _3/2 [/C (R, r, t) - /C (R, r, +oo)], (115) 

provided that R does not lie too close to dTZ. The limit -Ko(R> r ;+°°) is the inverse 
Fourier transform, g (R,,r), ofgo(Tl,k). Asymptotic formula (19T1) applies accordingly. 
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5. Conclusion 

We studied aspects of the time-dependent pair-excitation function, K , introduced by 
Wu [91 [10] for interacting Bosons at zero temperature. The particles are trapped by a 
sufficiently slowly varying external potential. By assuming that the wavefunction, $, of 
the macroscopic state satisfies a time-independent nonlinear Schrodinger equation, we 
determined an approximate solution for the time-dependent integrodifferential equation 
for Kq. Our analysis relied on the fact that, because of the slowly varying external 
potential, the space variables in K are separated into the fast variable r = x — y and 
the slow variable R = (x + y)/2 of the center-of-mass system. The ensuing solution for 
Kq is given in terms of the Fourier transform in r where R enters as a vector parameter. 

For zero pair excitation (K = 0) initially, R lying inside the trap and sufficiently 
large t, i.e. t ^> (lGnapo^oCR) 2 )' 1 where $o(R) is given by (1961) . we derived 
asymptotic formula fl9T|) . This formula describes how K approaches its steady-state 
values, go(R) r )- This result involves a convergent series containing Lommel's functions 
with arguments depending on the (scaled) variables f and r of (11141) . Denoting 
[167rap $o(R) 2 ] 1/2 by u(R), where u becomes the sound velocity when $o is unity |15j . 
the asymptotic analysis reveals that \K$ — go I takes small, 0(t~ A ) values for r <C 2ut 
but increases appreciably to become 0({r-t)~ r ) as r approaches 2ut. 

The condition r — n 2ut = 0((nt) 1/3 ), where n is any positive integer, plays a key 
role in our asymptotic analysis: it signifies that a stationary-phase point in one of the 
requisite integrals falls too close to the endpoint of integration. This effect is reminiscent 
of the coalescence of critical points and the onset of caustics in diffraction theory [19J, 
but the special functions involved are of course different. A physical interpretation of 
such coalescence in the present case is tempting but elusive. 

Our assumption of $ satisfying the nonlinear Schrodinger equation is a crucial one. 
Clearly, $ is modified when Kq acts back on it. This issue is not addressed by our 
analysis. Furthermore, the assumption of a time-independent, slowly varying external 
potential may have to be relaxed in order to account for a wider class of atomic traps. 
This consideration is left for future study. 
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Appendix A. Solution of differential equation for /C 



In this appendix we solve the ordinary differential equation (1241) . 
Consider the ordinary differential equation 

iz = z 2 + 2(a 2 + l)z + 1 , 



(A.l) 
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where z = z(t), a is a positive constant and the dot on top of z denotes differentiation 
with respect to t. By factorizing the right-hand side of (lA.lj) we have 



iz= (z + 1 + a 2 - oV2 + a 2 ){z + 1 + a 2 + aV2 + a 2 ) . (A.2) 
Separating variables, direct integration of the last equation yields the general solution 

i 

In 



2ay/2 



a 



2 



z + l + a 2 - ay/2 + 



z+ l + a 2 + ay/2 + 



Q 



which is in turn solved for z(t) to give 



t + C, (A.3) 



z(t) = ay/2 + a 2 - 1 - a 2 + 2a\/2 + a 2 ==- , (A.4) 

where C is an integration constant. 

The constant C is computed by applying the initial condition z(0) = zq. We find 



z + 1 + a 2 - aVl + 



a 



2 



C = ^ u ' ' ^ 7 Y ' . (A.5) 

1 + a 2 - aV2 + a 2 

Equation fl28|) is recovered by identifying a with (Snapo^ 1 / 2 k and z with /(k). 

Appendix B. Complex singularities of K for zero initial data 

In this appendix we discuss the analytic continuation of the right-hand side of fl32|) in 
the complex fc-plane. We focus on the simplified yet physically appealing case with 
p{k) = goik), which results from zero initial data, /(r) = in (1201) . 
In units where 167rapo — 1 the function of interest reads 

u(k) = Uk) + Mk) i! ^^-^ , (B.i) 



where u{k) = kyW+1 and 

g (k) = -2(k 2 + \ - u) , ? (A:)- 1 = -2(A; 2 + | + ^) . (B.2) 
After some straightforward algebra, we obtain the formula 

lik 1 + |J S[uj) + cos\uJt) 

where 

SM := . (B.4) 

Because S(u) is analytic and even in u>, the only possible singularities of U{k) are poles. 
Note that the branch points k = and ±i of u(k) are regular points of U(k). 

Next, we discuss the location of the poles of U (k). With the substitution k = sinh n , 
we find that the denominator in ( IB. II) vanishes at points 77 where 

|t sinh(2?7) — i 2r\ = mil, m = ±1, ±2, ... . (B.5) 

Let 7] = a + ib. The imaginary part of the last equation yields 

4a = t sin(26) cosh(2a) (t > 0) . (B.6) 
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It suffices to consider only the range —tt/2 < b < tt/2. It is readily concluded that it is 
impossible to have a > and —tt/2 < b < 0; or, a < and < b < tt/2. Thus, the poles 
must lie in the first and third quadrants of the complex fc-plane. A further examination 
of (IB. 51) shows that these poles are simple. 

It is of interest to describe the poles analytically for small and large values of \rmr\/t. 
For t ^> 1 and — n < 2b < ir, we find that flB.5j) is solved approximately by 

17171 . . 

77 ~ , \m\7T <C t, m =1,2... , 

t — 2% 'ii 

I / % \ 4 m 7r 
r? ~ - sg(m) I 1 + | — — I In , \m\n 3> t , (B-7) 



2 \ |m|7r / i 

where sg is the sign function, i.e. sg(x) = 1 if x > and sgn(x) = — 1 if x < 0. 
Evidently, the = sinh?7 corresponding to flB.71) lies in the first or third quadrant of 
the k plane. For fixed t, these poles approach the real axis. 

Appendix C. Analytic continuation formulas for S i(z) 

In this appendix, S i (z e~ l7T ) and S i(z e t2n ), which are involved in the asymptotics 
for A(r, t) in section |3T2| are expressed in terms of S Q i(z). 
The starting point is the relation [23j 

S i(z) = s 0>h (z) + |r(| - |) r(| + ±) [- sin(7r/6) J 1/3 (z) 
- cos(7r/6) Y 1/3 (z)] 
= - \ [tan(7r/6) J 1/3 (z) + Y 1/3 (z)} , (C.l) 

where T(z) := J °° dxx 2_1 e _a; is the Gamma function, J u (z) and 1^(2) are Bessel 
functions [23] . 

9^ 



s i(z) := - iF 2 (l;4/3, 5/3; -z*/4) , (C.2) 



and iF 2 is a hypergeometric series [24] . Evidently, we have 

s i (z e—) = - S0)| (2) = -s^_ (z e^) . (C.3) 

Useful analytic continuation formulas for the Bessel functions are [23] 

J 1/3 (ze-n=e~ t7r/3 J 1/3 (z) , J 1/3 (z e 12 *) = e i2 ^ 3 J 1/3 (z) , 
r 1/3 (^e- i7r ) = e^ 3 Y 1/3 (z) -1 J 1/3 (z) , 

Y 1/3 (z e^) = e"^ 3 y 1/s (z) + i J 1/3 (z) . (C.4) 
By combining (IC.1|) - (IC.4I) we obtain the desired analytic continuation formulas, 
Sn&e-**) = -S i(z) -^V3e-^ s H[%(z) , (C.5) 



- 2 1 :i( 

7T 



S Aze^) = S Uz) + I V3e-^ 3 H[%(z) , (C.6) 



where H^\z) is the Bessel function of the third kind 
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Equations (10.5j> and (10.61) hold for any z. Replacing z by % z we have 

S 0>1 _(ize-™) = - S i(iz) + V3K 1/3 (z) , (C.7) 

S i(ize i2 *) = S x(iz) - V3K 1/3 (z) , (C.8) 
where K v (z) is the modified Bessel function of the third kind |23j . 
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